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In a sample of 471 million BB events collected with the BABAR detector at the PEP-II e"'"e~ collider 
we study the rare decays B — )■ K^''^t^i~ , where i~ is either e'^e~ or fi'^fi' . We report results 
on partial branching fractions and isospin asymmetries in seven bins of di-lepton mass-squared. We 
further present CP and lepton-flavor asymmetries for di-lepton masses below and above the J/tp 
resonance. We find no evidence for CP or lepton-flavor violation. The partial branching fractions 
and isospin asymmetries are consistent with the Standard Model predictions and with results from 
other experiments. 

PACS numbers: 13.20.He 



I. INTRODUCTION 

The decays B K^*H+l- arise from flavor-changing 
neutral-current processes that are forbidden at tree level 
in the Standard Model (SM). The lowest-order SM pro- 
cesses contributing to these decays are the photon pen- 
guin, the Z penguin and the W~^W~ box diagrams shown 
in Fig. [1] Their amplitudes are expressed in terms of 
hadronic form factors and perturbatively-calculable ef- 
fective Wilson coefficients, C^^ , and CJq , which 
represent the electromagnetic penguin diagram, and the 
vector part and the axial- vector part of the linear combi- 
nation of the Z penguin and W^W~ box diagrams, re- 
spectively Jj. In next-to- next-to-leading order (NNLO) 
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at a renormalization scale /i = 4.8 GeV, the effective 
Wilson coefficients are Cf = -0.304, Cf = 4.211, and 
Cfff = -4.103 i. 

Non-SM physics may add new penguin and box dia- 
grams, which can contribute at the same order as the SM 
diagrams 0-13. Examples of new physics loop processes 
are depicted in Fig. [2] These contributions might modify 
the Wilson coefficients from their SM expectations iSt-tSj] . 
In addition, new contributions from scalar, pseudoscalar, 
and tensor currents may arise that can modify, in partic- 
ular, the lepton-flavor ratios [ol. [lol|. 




q q q q 



FIG. 1: Lowest-order Feynman diagrams for b se+r. 
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(a) W 



(b) %- 



(C) S:t 



metry is defined as 



b t,c,u s b ~t,c,u s b b,s,d s 

FIG. 2: Examples of new physics loop contributions to 
b — !> s£'^£^: (a) charged Higgs {H~); (b) squark {i,c,u) and 
chargino (c) squark {b,s,d) and gluino (gt)/neutralino 

(x°). 



(.) 
CP 



Ki*)(+t 



B{B K^*\+l-) + B{B Ki*)£+i-) 



(1) 



and is expected to be 0(10^"^) in the SM. However A^p ' 
may receive a significant enhancement from new physics 
contributions at the electro- weak scale [l^ . 

For s > 0.1 GeV^/c^, the ratio of rates to dimuon and 
dielectron final states is defined as 



II. OBSERVABLES 

We report herein results on exclusive partial branching 
fractions and isospin asymmetries in six bins of s = m^^, 
defined in Table HI We further present results in the s bin 
So — 1.0 — 6.0 GeV^/c^ chosen for calculations inspired 
by soft-coUinear effective theory (SCET) [ll|. In addi- 
tion, we report on direct CP asymmetries and the ratio 
of rates to dimuon and dielectron final states in the low s 
and high s regions separated by the J/V' resonance. We 
remove regions of the long-distance contributions around 
the J/i) and V'(2S') resonances. New BABMi results on 
angular observables using the same dataset and similar 
event selection will be reported shortly. 

The B Ki+£^ and B !<*£+£- total branch- 
ing fractions are predicted to be (0.35 ± 0.12) x 10'^ 
and (1.19 ± 0.39) x lO''^ (for s > 0.1 GeVVc"*), re- 
spectively [Sj. The ~ 30% uncertainties are due to a 
lack of knowledge about the form factors that model the 
hadronic effects in the B K and B K* transi- 
tions. Thus, measurements of decay rates to exclusive 
final states are less suited to searches for new physics 
than rate asymmetries, where many theory uncertainties 
cancel. 

For charged B decays and neutral B decays flavor- 
tagged through K* ->• JsT+tt" [i1|, the direct CP asym- 



n 



^ B{B ^ K'^*^H+^i-) 
" B{B ^ K(*)e+e-) 



(2) 



In the SM, TZj^m is expected to be unity to within 



few percent [lj| for dilepton invariant masses above the 
dimuon kinematic threshold. In two-Higgs-doublet mod- 
els, including supersymmetry, these ratios are sensitive 
to the presence of a neutral Higgs boson. When the ratio 
of neutral Higgs field vacuum expectation values tan (3 is 
large, TZ^i*) might be increased by up to 10% [l3|. 
The CP-averaged isospin asymmetry is defined as 



B{B'^^K(*'>0£+£-)+r-rB{B+^K(*'> + £+£-) ' 



(3) 

where = TBo/Tp+= 1/(1.071±0.009) is the ratio of 5° 
and 5+ lifetimes [la. Af has a SM expectation of -1-6% 
to -fl3% as s — )■ This is consistent with the mea- 
sured asymmetry 3 ±3% in _B — )■ K*j [l^. A calculation 
of the predicted K*'^ and K*''' rates inte grat ed over the 
low s region yields Af' = -0.005 ± 0.020 ]l3, [11. In the 
high s region, we may expect contributions from charmo- 
nium states as an additional source of isospin asymmetry. 
However the measured asymmetries in the J/ipK^*^ and 
V'(2S')X(*) modes are all below 5% [H. 



III. BABAR EXPERIMENT AND DATA 
SAMPLE 



TABLE I: The definition of seven s bins used in the anal- 
ysis. Here ms and rnj^(,) are the invariant masses of 
B and K^*\ respectively. The low s region is given by 
0.10 < s < 8.12 GeV^/c*, while the high s region is given 
by s > 10.11 GeV^c^ 





s bin 


s min 
( GeV^c") 


s max 
(GeVVc-*) 


Low 


Sl 


0.10 


2.00 




S2 


2.00 


4.30 




S3 


4.30 


8.12 


High 


S4 


10.11 


12.89 




S5 


14.21 


16.00 




Sd 


16.00 


{niB - mj^(,)f 




SO 


1.00 


6.00 



We use a data sample of 471 million BB pairs col- 
lected at the T{AS) resonance with the BaBAR detec- 
tor [l9| at the PEP-II asymmetric-energy e'^e" collider 
at the SLAG National Accelerator Laboratory. Charged 
particle tracking is provided by a five-layer silicon vertex 
tracker and a 40-layer drift chamber in a 1.5 T solenoidal 
magnetic field. We identify electrons with a CsI(Tl) 
electromagnetic calorimeter, and muons using an instru- 
mented magnetic flux return. Electron and muon candi- 
dates are required to have momenta p > 0.3 GeV/c in the 
laboratory frame. We combine up to three photons with 
electrons when they are consistent with bremsstrahlung, 
and do not use electrons that are associated with photon 
conversions to low-mass e'^e~ pairs. We identify charged 
kaons using a detector of internally reflected Cherenkov 
light, as well as dE/dx information from the drift cham- 
ber. Charged tracks other than identified e, ^ and K 
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candidates are treated as pions. Neutral K'^ — tt+tt" 
candidates are required to have an invariant mass con- 
sistent with the nominal K'^ mass, and a flight distance 
from the e"'"e~ interaction point that is more than three 
times its uncertainty. 



IV. EVENT SELECTION 

We reconstruct B — >■ K^*H~^£~ signal events in the 
following eight final states: 



• 5° ^ K2e+e'. 



e ' e 



• B° ^ K*"{-^ K+7r-)iJ,+fj.-, 

• B+ ^ K*+{-^ K2:-K+)e+e-, 

• B° ^ K*"{-^ K+n-)e+e-. 

We reconstruct K'^ candidates in the tt+tt" final state. 
We also study the K'^*^h^ix^ final states, where h is 
a charged track with no particle identification require- 
ment applied, to characterize backgrounds from hadrons 
misidentified as muons. We use a K*e^^^ sample to 
model the combinatorial background from two random 
leptons. In each mode, we utilize the kinematic vari- 
ables mEs = and AE ^ E*g - Ecm/2, 
where p*g and Eg are the B momentum and energy in 
the T{AS) center-of-mass (CM) frame, and Eqm is the 
total CM energy. 

For masses rriES > 5.2 GeV/c^ we perform one- 
dimensional fits of the mEs distribution for Ki^£~ 
modes. For K*£^£~ modes, we include in addition the 
Ktt mass region 0.72 < rriKn < 1-10 GeV/c^ in the fit. 
We use the sideband 5.20 < jties < 5.27 GeV/c^ to char- 
acterize combinatorial background shapes and normal- 
izations. For both the e~^e~ and ^.^ modes, we veto 
the J/V'(2.85 < m« < 3.18 GeV/c^) and V(25)(3.59 < 
mu < 3.77 GeV/c^) mass regions. The vetoed events 
provide high-statistics control samples that we use to val- 
idate the fit methodology. 

The main backgrounds arise from random combina- 
tions of leptons from semileptonic B and D decays. 
These combinatorial backgrounds from either BB events 
(referred to as "BS backgrounds") or continuum qq 
events (e+e^ qq, q — u, d, s, c, referred to as "gg 
backgrounds" ) are suppressed using bagged decision trees 
(BDTs) [205. We train eight separate BDTs as follows: 

• Suppression of BB backgrounds for e+e^ modes in 
the low s region; 



• Suppression of BB backgrounds for e+e modes in 
the high s region; 

• Suppression of BB backgrounds for ^jl^ ^jl^ modes 
in the low s region; 

• Suppression of BB backgrounds for /u^/i^ modes 
in the high s region; 

• Suppression of qq backgrounds for e~ modes in 
the low s region; 

• Suppression of qq backgrounds for e^e~ modes in 
the high s region; 

• Suppression of qq backgrounds for ^^fj,^ modes in 
the low s region; 

• Suppression of qq backgrounds for /i^/i^ modes in 
the high s region. 

The BDT input parameters include the following observ- 
ables: 

• AE of the B candidate; 



• The ratio of Fox- Wolfram moments R2 [2l| and the 
ratio of the second-to-zeroth angular moments of 
the energy flow L2/L0 113 1 both event shape pa- 
rameters calculated using charged and neutral par- 
ticles in the CM frame; 

• The mass and AE of the other B meson in the event 
(referred to as the "rest of the event" ) computed in 
the laboratory frame by summing the momenta and 
energies of all charged particles and photons that 
are not used to reconstruct the signal candidate; 

• The magnitude of the total transverse momentum 
of the event in the laboratory frame; 

• The probabilities that the B candidate and the 
dilepton candidate, respectively, originate from a 
single point in space; 

• The cosine values of four angles: the angle between 
the B candidate momentum and the beam axis, the 
angle between the event thrust axis and the beam 
axis, the angle between the thrust axis of the rest of 
the event and the beam axis, and the angle between 
the event thrust axis and the thrust axis of the rest 
of the event, all defined in the CM frame. 

Figure [3] shows the output distributions of the BDTs for 
Monte Carlo (MC) simulated signal and combinatorial 
background for the e+e^ sample below the J/ip reso- 
nance. The distributions are histograms normalized to 
unit area. The selections on BDT outputs are further 
optimized to maximize the statistical significance of the 
signal events, as shown later. 

Another source of background arises from B D{-^ 
K''*^ 71)77 decays if both pions are misidentified as lep- 
tons. Determined from data control samples with high 
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FIG. 3: The (a) BB and (b) qq e+e" BDT outputs for 
simulated events in the low s region. Shown are the distri- 
butions for BB background (red dashed line), qq background 
(red dotted line), and signal (blue solid line) event samples, 
normalized to unit area. 



purity pL^ , the misidentification rates for muons and elec- 
trons are ^ 3% and 0.1% per candidate, respectively. 
Thus, this background is only significant for fi- 
nal states. We veto these events by requiring the in- 
variant mass of the K'^*^tt system to be outside the range 
1.84—1.90 GeV/c^ after assigning the pion mass hypothe- 
sis to the muon candidates. Any remaining residual back- 
grounds from this type of contribution are parameterized 
using control samples obtained from data. 

After applying all selection criteria about 85% of signal 
events contain more than one B candidate. These can- 
didates differ typically in one charged or neutral hadron. 
The average number of candidates per signal event is 
about six. To choose the best candidate, we define the 
ratio 



A = 



V 



BB 
sig 



-pit 

sig 



V 



BB 



■pqq 

sig 



■pBB 
I bkg 



■pqq ■ 

' bkg 



(4) 



where Vsig and Vbkg are probabilities calculated from 
the corresponding BB and qq BDT output distributions 



for signal and background, respectively. We select the 
candidate with the largest A as the best candidate. The 
probability for a correctly-reconstructed signal event to 
be selected as the best candidate is mode-dependent and 
varies between about 80% and 95% for s bins below the 
J/t/j mass, while for s bins above the ^^{23) mass it varies 
between about 60% and 90%. 



V. SELECTION OPTIMIZATION 

To optimize the AE selection, we simultaneously vary 
the upper and lower bounds of the AE interval to find the 
values that maximize the ratio S/y/S + B in the signal 
region (mEs > 5.27 GeV/c^, and for K* modes in addi- 
tion 0.78 < rriK-nr < 0.97 GeV/c^), where S and B are 
the expected numbers [l5l | of signal and combinatorial 
background events, respectively. We perform separate 
optimizations for dilepton masses below and above the 
J/f/) mass. For some modes, the optimization tends to 
select very narrow intervals, which leads to small signal 
efficiency. To prevent this, we require the magnitudes 
of the AE upper and lower bounds to be 0.04 GeV or 
larger. (Note that the lower bound is always negative 
and the upper bound always positive.) 

We also optimize the lower bounds on the BDT BB 
and qq intervals (the upper bounds on these intervals are 
always 1.0). We perform fits to extract signal yields using 
the fit model described in Sec. IVII For each mode, the 
lower bound on the BDT interval is optimized by max- 
imizing the expected signal significance defined as the 
fitted signal yield divided by its associated uncertainty. 
We determine these from 500 pseudo-experiments using 
branching fraction averages [15j . The optimized BDT 
lower bounds are listed in Tables |TT] and IIIII for K£^£^ 
and K*i^£~, respectively. Figure |4] shows the expected 
experimental significance in the BB BDT versus the qq 
BDT plane for B^ — )■ K^'k~ fX^ ii~ in bin S2- The signal 
selection efficiency and the cross-feed fraction (defined in 
Sec. IVII) in each mode and s bin after the final event selec- 
tion are also listed in Tables Inland IIIII The selection effi- 
ciencies determined in simulations vary from 11.4±0.2% 
for i^T^TT+e+e" in sg to 33.3 ± 0.3% for K^fi^fj.^ in ss, 
where the uncertainties are statistical. 



VI. FIT METHODOLOGY 

We perform one-dimensional fits in tties for Kt^l~ 
modes and two-dimensional fits in tties and rriK-n for 
K*t^i~ modes to extract the signal yields. The probabil- 
ity density function (PDF) for signal toes is parametrized 
by a Gaussian function with mean and width fixed to 
values obtained from fits to the vetoed J /ip events in the 
data control samples. For rriK-K, the PDF is a relativis- 
tic Breit-Wigner line shape [23|. True signal events are 
those where all generator-level final-state daughter parti- 
cles are correctly reconstructed and are selected to form 
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FIG. 4: Expected statistical significance of the number of fit- 
ted signal events as a function of BDT interval lower bounds 
for B° K+Tv-fi+^i- in bin S2- The star marks the opti- 
mized pair of lower bounds. 



TABLE II: Optimized lower bounds on the BDT intervals, 
signal reconstruction efficiency, and cross-feed fraction, by 
K£^£~ mode and s bin. The uncertainties are statistical only. 





Mode 


s bin 


BB 




Efficiency 


Cross- feed 








BDT 


BDT 


[%] 


fraction [%] 






Sl 


0.20 


0.80 


19.9 ± 0.2 


8.9 ±0.3 






S2 


0.70 


0.85 


22.2 ± 0.2 


8.6 ±0.2 






S3 


0.20 


0.85 


25.2 ±0.1 


8.9 ±0.2 






S4 


0.70 


0.70 


24.3 ±0.2 


9.4 ±0.2 






S5 


0.70 


0.80 


22.2 ±0.2 


12.0 ±0.5 






S6 


0.75 


0.80 


16.6 ±0.1 


21.7 ±0.7 






So 


0.50 


0.85 


22.7 ±0.1 


8.8 ±0.1 


B+ 


-> K+fi+fi- 


Sl 


0.30 


0.85 


21.3 ± 0.2 


0.3 ±0.0 






S2 


0.15 


0.85 


27.0 ± 0.2 


0.3 ±0.0 






S3 


0.15 


0.85 


30.9 ±0.1 


0.3 ±0.0 






S4 


0.80 


0.85 


31.0 ±0.2 


0.4 ±0.0 






S5 


0.65 


0.85 


33.3 ±0.3 


2.1 ±0.1 






S6 


0.05 


0.85 


30.5 ±0.2 


10.4 ±0.2 






So 


0.05 


0.85 


13.6 ±0.1 


0.3 ±0.0 






Sl 


0.25 


0.80 


22.1 ± 0.2 


8.3 ±0.3 






S2 


0.25 


0.80 


25.2 ± 0.2 


9.4 ±0.3 






S3 


0.65 


0.80 


24.3 ±0.1 


9.4 ±0.2 






S4 


0.50 


0.85 


24.1 ±0.2 


10.9 ±0.4 






S5 


0.05 


0.65 


23.0 ±0.2 


18.5 ±0.9 






S6 


0.25 


0.70 


16.5 ±0.1 


35.0 ±1.1 






So 


0.85 


0.85 


21.3 ±0.1 


9.2 ±0.2 


B+ 




Sl 


0.35 


0.85 


22.8 ±0.2 


0.4 ±0.1 






S2 


0.10 


0.85 


28.8 ±0.2 


0.4 ±0.0 






S3 


0.10 


0.85 


30.8 ±0.1 


0.5 ±0.0 






S4 


0.30 


0.80 


32.7 ±0.2 


1.1 ±0.1 






S5 


0.25 


0.80 


31.7 ±0.3 


4.3 ±0.2 






S6 


0.50 


0.85 


25.1 ±0.2 


12.0 ±0.3 






So 


0.40 


0.85 


29.6 ±0.1 


0.5 ±0.0 



TABLE III: Optimized lower bounds on the BDT intervals, 
signal reconstruction efficiency, and cross-feed fraction, by 
K*t^ l~ mode and s bin. The uncertainties are statistical 
only. 



Mode 






s bin 


BB 


qq 


Efficiency 


Cross-feed 










BDT BDT 


[%] 


fraction [%] 


B+ ^ 




MM 


Sl 


0.55 


0.85 


13.6 ±0.1 


14.0 ± 0.5 








S2 


0.80 


0.85 


14.6 ±0.2 


19.2 ± 0.7 








S3 


0.85 


0.80 


14.9 ±0.1 


20.7 ±0.5 








S4 


0.85 


0.85 


14.7 ±0.1 


28.0 ±0.7 








S5 


0.15 


0.85 


16.4 ±0.2 


59.3 ± 1.3 








S6 


0.10 


0.85 


14.3 ±0.1 


110.8 ± 1.9 








So 


0.80 


0.85 


14.5 ±0.1 


18.9 ±0.5 




K+-K- 




Sl 


0.80 


0.85 


16.2 ±0.1 


4.9 ±0.2 








S2 


0.80 


0.85 


19.6 ±0.2 


7.8 ±0.3 








S3 


0.75 


0.85 


21.3 ±0.1 


10.1 ±0.2 








S4 


0.85 


0.85 


20.9 ±0.1 


13.8 ± 0.3 








S5 


0.75 


0.85 


22.8 ±0.2 


31.7 ±0.6 








S6 


0.80 


0.80 


19.5 ±0.2 


61.0 ±0.9 








So 


0.60 


0.85 


20.4 ±0.1 


8.9 ±0.2 


B+^ 


kW 


- + - 
e e 


Sl 


0.45 


0.70 


16.6 ±0.2 


17.8 ± 0.6 








S2 


0.85 


0.85 


13.7 ±0.2 


20.7 ±0.8 








S3 


0.55 


0.85 


16.0 ±0.1 


27.5 ± 0.7 








S4 


0.40 


0.85 


15.4 ±0.1 


41.6 ±0.9 








S5 


0.80 


0.45 


13.1 ±0.2 


68.6 ± 1.8 








S6 


0.60 


0.85 


11.4 ±0.2 


133.4 ± 2.9 








So 


0.70 


0.85 


16.0 ±0.1 


23.1 ±0.5 


B" 


K+-K- 




Sl 


0.80 


0.85 


16.5 ±0.2 


6.8 ±0.2 








S2 


0.85 


0.85 


18.6 ±0.2 


10.9 ± 0.3 








S3 


0.80 


0.80 


18.5 ±0.1 


11.2 ±0.3 








S4 


0.55 


0.65 


21.9 ±0.2 


25.6 ±0.4 








S5 


0.75 


0.80 


19.0 ±0.2 


50.4 ±0.9 








S6 


0.05 


0.80 


15.1 ±0.2 


110.9 ± 1.8 








So 


0.80 


0.85 


19.7 ±0.1 


10.8 ±0.2 



a B candidate. 

For the combinatorial background, the mEs PDF 
is modeled with a kinematic threshold function whose 



shape is a free parameter in the fits [2Jj, while the rriK-K 



PDF shape is characterized with the K*e fi^ sample 



mentioned in Sec. IIVI We parameterize the combina- 
torial rriKTv distributions with non-parametric Gaussian 
kernel density estimator shapes [2q| (referred to as the 
"KEYS PDFs") drawn from the i^*e=^/z=F sample in the 
full Toes fit region. Since the correlation between mpcTr 
and AE is weak, we accept all K*e^fi^ events within 
\AE\ < 0.3 GeV, rather than imposing a stringent AE 
selection, in order to enhance sample sizes. 

Signal cross-feed consists of mis-reconstructed signal 
events, in which typically a low-momentum tt^ or tt'^ is 
swapped, added, or removed in the B candidate recon- 
struction. We distinguish among different categories of 
cross-feed: "self-cross-feed" is when a particle is swapped 
within one mode, "feed-across" is when a particle is 
swapped between two signal modes with the same final- 
state multiplicity, and "feed-up (down)" is when a parti- 
cle is added (removed) from a lower (higher) multiplicity 
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b s£~^£~ mode. We use both exclusive and inclusive 
b — > s£^£^ MC samples to evaluate the contributions 
of the different categories. The cross-feed niES distri- 
bution is typically broadened compared to correctly re- 
constructed signal decays. We combine the cross-feed 
contributions from all sources into a single fit compo- 
nent that is modeled as a sum of weighted histograms 
with a single overall normalization, which is allowed to 
scale as a fixed fraction of the observed correctly recon- 
structed signal yield. This fixed fraction is presented as 
the "cross-feed fraction" in Tables and IIIIl The mod- 
eling of cross-feed contributions is validated using fits to 
the vetoed J/ipK^*^ and 'iIj{2S)K'^*'^ events, in which the 
cross-feed contributions are relatively large compared to 
all other backgrounds. 

Exclusive B hadronic decays may be mis-reconstructed 
as , since hadrons can be misidentified as 

muons. Following a procedure similar to that described 
in Ref . 26] , we determine this background by selecting a 
sample of K'^*^ ^^h^ events, in which the muon is iden- 
tified as a muon and the hadron is inconsistent with 
an electron. Requiring identified kaons and pions, we 
select subsamples oi K'^*^t:+it^ ,K''*^K+tt~' ,K^*'>'n:+K'' , 
and K^*^K+K~. We obtain weights from data control 
samples where a charged particle's species can be identi- 
fied with high precision and accuracy without using par- 
ticle identification information. The weights are then ap- 
plied to this dataset to characterize the contribution ex- 
pected in our fits due to misidentified muon candidates. 
We characterize the misidcntification backgrounds using 
the KEYS PDFs, with normalizations obtained by con- 
struction directly from the weighted data. 

Some charmonium events may escape the charmonium 
vetoes and appear in our fit region. Typically, this oc- 
curs when electrons radiate a photon or a muon candi- 
date is a misidentified hadron and the missing energy is 
accounted for by a low-energy ir^ or tt". The largest 
background contributions from this source are expected 
in the K* ji^ ji'^ and _ftr*e+e^ channels. We model this 
background using the charmonium MC samples and de- 
termine the leakage into s bins on either side of the J / 
and ip{2S) resonances. We see a notable charmonium 
contribution (about five events) for B^ — )■ K~^Tr~ fi~ 
in bin S3. This leakage is typically caused by a swap be- 
tween the /i"*" and tt"*" in a single B — >■ J/tp (— >■ fX^ ijl~)Ktt^ 
candidate, where both the fX^ and tt"*" are misidentified. 

Hadronic peaking background from B — K*Tr^ and 
B — > K*ri in which the 7r° or 77 decays via Dalitz pairs 
shows a small peaking component in mes in bin si. Due 
to the requirement s > 0.1 GeV^/c^, contributions of 7 
conversions from B — > K*"f events beyond the photon 
pole region are found to be negligible. 



Fit Model for Rate Asymmetries 

Using the PDFs described above, we perform simul- 
taneous fits across different K'-*k+£- modes, bmce 



efficiency-corrected signal yields are shared across vari- 
ous decay modes, we can extract rate asymmetries di- 
rectly from the fits. The fitted signal yields in B^ modes 
are corrected by the lifetime ratio Tga /tb+- We also cor- 
rect the signal yields for B{K* — ^ Ktt) in K* modes and 
B{K^ — > TT+TT^) in the modes with a K^. In the fits for 
Acp, we share the efficiency-corrected signal yield Nb as 
a floating variable for B (qb, q = u,d) events across differ- 
ent flavor-tagging iir(*'£+£~ modes by assuming lepton- 
flavor and isospin symmetries. The efficiency-corrected 
signal yield Ng for B {qb) events is then deflned by 
Ng = Nb • (1 + -^CP))/(1 ~ -^cp) and is also shared 
across corresponding modes. For the lepton-flavor ra- 
tios TZji(,) , we share the efliciency-corrected signal yield 
Nee as a floating variable for the two B — > Ke'^e~ or 
B — > K*e'^e~ modes by assuming isospin symmetry. 
The efficiency-corrected signal yield A^^^ shared across 
the corresponding B — > modes is then de- 

fined by 7V^^ = Nee ■ Rk'*)- Foi' the isospin asymme- 
try Af'' ' , we share the efficiency-corrected signal yield 
Nb+ as a floating variable for the two — K'^£^£~ or 
B^ — > K*^l'^£~ modes by assuming lepton-flavor sym- 
metry. The efficiency-corrected signal yield iV^o shared 
across the corresponding i?" — > K*^£'^£'^ modes is then 
deflned by Nb^ = Nb+ ■ (1 + Af^*^)/{1 - ylf 



VII. FIT VALIDATION 

We validate the fit methodology with charmonium con- 
trol samples obtained from the dilepton mass regions 
around the J/tp and 'il'{2S) resonances that are vetoed in 
the B K'-*'>£+i~ analysis. We measure the J/ipK'-*'> 
and ^p{2S)K^*^ branching fractions in each final state 
with the optimized BDT selections in bins S3 and S4, re- 
spectively. Our measurements agree well with the world 
averages [lB| for all final states. Typical deviations, 
based on statistical uncertainties only, are less than one 
standard deviation (a). The largest deviation, in the 
Tr~ fj,'^ IJ,~ mode, is 1.7a. For J/tpK^*^ modes, the sta- 
tistical uncertainties are considerably smaller than those 
of the world averages. We fioat the Gaussian means and 
widths of the signal PDFs in the fits for the J/tJjK^*^ 
modes. The associated uncertainties obtained from the 
fits arc then used as a source of systematic variation for 
the signal PDFs. The typical signal width in rngs is 
2.5 MeV/c2. 

We further validate our fitting procedure by applying 
it to charmonium events to extract the rate asymme- 
tries. The measured CP asymmetries Acp, lepton-flavor 
ratios Ti-x(*) and isospin asymmetries Ai are in good 
agreement with Standard Model expectations or world 
averages for Ai . 

We also test the methodology with flts to ensem- 
bles of datasets where signal and background events are 
generated from appropriately normalized PDFs ("pure 
pseudo-experiments"). We perform flts to these pseudo- 
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experiments in each mode and s bin using the full fit 
model described previously. For ensembles of 1000 pure 
pseudo-experiments, the pull distributions for the sig- 
nal yields show negligible biases. We further fit ensem- 
bles of pseudo-experiments in which the signal events are 
drawn from properly normalized exclusive MC samples 
("embedded pseudo-experiments"). The pull distribu- 
tions also show the expected performance. 

We perform fits to ensembles of pure pseudo- 
experiments in order to estimate the statistical sensi- 
tivity of, and biases related to, the various rate asym- 
metry measurements. The pull distributions for Acp 
and TZj^M for the low and high s regions show mini- 
mal biases. For Ai, we test a series of Ai input val- 
ues (—0.6, —0.3, 0.0, 0.3, 0.6) in each s bin using pure 
pseudo-experiments to ensure we obtain unbiased fits un- 
der different assumptions of isospin asymmetry. The Af 
pulls are generally well-behaved. In the worst case, the 
test fits for A^ are slightly biased due to very low signal 
yield expectations in the K^i'^l'^ final states. 

VIII. SYSTEMATIC UNCERTAINTIES 

Since some systematic uncertainties largely cancel in 
ratios, it is useful to separate the discussion of systematic 
uncertainties on partial branching fractions from that on 
rate asymmetries. 

A. Branching Fraction Uncertainties 

Systematic uncertainties for branching fractions arise 
from multiplicative systematic uncertainties involving 
the determination of the signal efficiency, and from addi- 
tive systematic uncertainties arising from the extraction 
of signal yields in the data fits. The multiplicative sys- 
tematic errors include contributions from the 

• Number of BB pairs: This uncertainty is 0.6%. 

• Tracking efficiency for charged particles: We as- 
sign a correlated uncertainty of 0.3% for each lep- 
ton, and 0.4% for each charged hadron including 
daughter pions from decay |27j . 

• Charged particle identification (PID) efficiencies: 
We employ a data-driven method to correct PID 
efficiencies in simulated events. We estimate the 
systematic uncertainties from the change in signal 
efficiency for simulated J/ipK^*^ events after turn- 
ing off the PID corrections. The systematic un- 
certainties are mode dependent and vary between 
0.3% and 1.6%. 

• ifg identification efficiency: This is determined as 
a function of fiight distance after applying if^ ef- 
ficiency corrections. An uncertainty of 0.9% is ob- 
tained by varying the K^^ selection algorithm. 



• Event selection efficiency: We measure the effi- 
ciency of the BDT selection in charmonium data 
control samples and compare with results obtained 
for exclusive charmonium samples from simulation. 
We take the magnitude of the deviation for any par- 
ticular final state and s bin as the uncertainty asso- 
ciated with the BDT lower bounds. If the data and 
simulation are consistent within the uncertainty, 
we then take the uncertainty as the systematic un- 
certainty. The systematic uncertainty is found to 
vary between 0.3% and 9.1% depending on both the 
mode and the s bin. Due to a strong correlation be- 
tween the A£' and BDT outputs, uncertainties due 
to Ai? are fully accounted for by this procedure. 

• Monte Carlo sample size: We find the uncertainty 
related to the finite size of the MC sample to be of 
the order of 1% or less for all modes. 

The additive systematic uncertainties involve contri- 
butions from the 

• Signal PDF shapes: We characterize them by vary- 
ing the PDF shape parameters (signal mean, signal 
width, and combinatorial background shape and 
normalization) by the statistical uncertainties ob- 
tained in the fits to the J/^ data control samples 
for Toes a-nd signal MC events for mK-n- 

• Hadronic backgrounds: We characterize them by 
varying both the normalization by the associated 
statistical uncertainties and by performing fits with 
different choices of smoothing parameters for the 
KEYS PDF shapes. 

• Peaking backgrounds from charmonium events and 
vr^/r/ Dalitz decays: We vary the normalization for 
these contributions by ±25%. 

• Modeling of to^c^ line shapes of the combinatorial 
background: We characterize the uncertainties by 
analyzing data samples selected from the toes < 
5.27 GeV/c^ sideband, and simulated events. 

Table IIVI summarizes all sources of systematic uncer- 
tainties considered in the total branching fraction mea- 
surements for individual modes. The total systematic 
uncertainty for the branching fractions is obtained by 
summing in quadrature the above-described uncertain- 
ties from different categories. 

B. Systematic uncertainties for the rate 
asymmetries 

For Acp^ a large portion of the uncertainties associated 
with the signal efficiency cancel. We find that the only 
efficiency- related term discussed in Sec. IVIIIAl that is not 
negligible for Acp is the one associated with the PID se- 
lection. Amongst the efficiency-related systematics, we 
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TABLE IV: Individual systematic uncertainties [%] for measurements of the total branching fractions in K^*^l^£ decays. 



Mode 






K\le+e- 


K+e+e- K'^TT- 






e+e K+n 


e+e 


BB counting 


±0.6 


±0.6 


±0.6 


±0.6 


±0.6 


±0.6 


±0.6 


±0.6 


Tracking 


±1.4 


±1.0 


±1.4 


±1.0 


±1.8 


±1.4 


±1.8 


±1.4 


PID 


±1.6 


±0.3 


±0.7 


±0.4 


±1.5 


±0.3 


±0.5 


±1.2 


RO id 


±0.9 




±0.9 




±0.9 




±0.9 




BDT selections 


±2.2 


±1.7 


±4.7 


±1.5 


±8.3 


±2.5 


±9.1 


±2.7 


MC sample size 


±0.3 


±0.3 


±0.3 


±0.3 


±0.4 


±0.3 


±0.4 


±0.4 


Sig. Shape 


±0.5 


±0.4 


±1.5 


±0.4 


±1.5 


±0.7 


±1.5 


±0.7 


Hadronic 


±3.3 


±5.8 






±2.3 


±1.6 






Peaking 


±0.3 


±0.8 


±1.2 


±0.8 


±0.7 


±1.7 


±0.8 


±1.2 


Comb. m^jT shape 










±1.2 


±0.6 


±0.6 


±1.6 


Total 


±4.7 


±6.3 


±5.4 


±2.2 


±9.3 


±3.9 


±9.5 


±4.0 



therefore only consider this term. We also consider the 
additive systematic uncertainties listed in Sec. IVIII Al 
Our measured Acp central values for J/^ K and J/tp K* 
are both well below 1% and show minimal detector effi- 
ciency effects. Potential, additional Acp systematic ef- 
fects from the assumptions of lepton-flavor and isospin 
symmetry are tested by removing these assumptions. 

The systematic uncertainties for the lepton-flavor ra- 
tios TZj^{,) are calculated by summing in quadrature the 
systematic errors in the muon and electron modes. Com- 
mon systematic effects, such as tracking, K'^ efficiency, 
and BB counting, yield negligible uncertainties in the 
ratios. Potential, additional Tlj^{*) systematic effects are 
tested by removing the assumption of isospin symmetry. 

For the systematic uncertainties of Ai, we sum in 
quadrature the systematic errors in charged and neutral 
B modes. Common systematic effects, which include 
BB counting and a large portion of the uncertainties 
associated with PID and tracking efficiencies, are neg- 
ligible. Again, additional tests on Aj systematics are 
performed by relaxing the assumption of lepton-flavor 
symmetry. Furthermore, as the cross-feed fractions in 
Tables HIl and IIIII are estimated under the assumption 
of isospin symmetry, we test this systematic effect using 
cross-feed fractions estimated with different Aj input val- 
ues. 

Our checks on symmetry assumptions described above 
for Acp, 'R-Ki-) and Ai generally show deviations from 
the original measured values below 20% of the associ- 
ated statistical uncertainties, and so we do not assign 
additional uncertainties. 



IX. RESULTS 

We perform fits for each K^*^t^l'' final state in each s 
bin listed in Tables |lT] and IIIII to obtain signal and back- 
ground yields, A^sig and A'^bkg, respectively. We model the 
different background components by the PDFs described 
in Sec. I VII We allow the shape parameter of the toes 
kinematic threshold function of the combinatorial back- 
ground to float in the fits. For the signal, we use a flxed 
Gaussian shape unique to each final state, as described 
previously. We leave the shapes of the other background 



PDFs fixed. For the peaking background, we fix the ab- 
solute normalization. For the cross-feed, we fix the nor- 
malization relative to the signal yields. 

Figure [5] shows as an example the toes distribution 
for the combined K£~^£~ modes in bin 54, while Fig. [S] 
shows the toes and mxTr mass spectra for the combined 
K*£'^£~ modes in bin si. The cross-feed contributions 
and the peaking backgrounds are negligible for this fit. 
The combinatorial background dominates and for /i^/i^ 
modes misidentified hadrons are the second largest back- 
ground. From the yields in each s bin we determine the 
partial branching fractions summarized in Table [V] Fig- 
ure [7] shows our results for the partial branching frac- 
tions of the K£^£^ and K*£^£^ modes in comparison to 
results from the Belle and CDF Collaborations 0, [29j 
and to the prediction of the Ali et al. model [HI ■ Our re- 
sults are seen to agree with those of Belle and CDF. Our 
results are also in agreement with the most recent par- 
tial branchingfraction measurements of B'^ — > if 
from LHCb ^J. 

The total branching fractions are measured to be 

BiB-^K£+£-) = (4.7±0.6±0.2) X 10"^ 
B{B ^ K*£+£-) = (10.21};;^ ± 0.5) X 10"^. 

Here, the first uncertainties are statistical, and the sec- 
ond are systematic. The total branching fractions are 
shown in Fig. [8] in comparison to measurements from 
Belle [m and CDF [29| and predictions from Ali et al. 3 
and Zhong et al. 

To measure direct Acp, we fit the B and B samples in 
the two K~^£~^£~ modes and four K*£~^£~ modes listed 
in Sec. lIVI We perform the measurements in the full s re- 
gion, as well as in the low s and high s regions separately. 
The B and B data sets share the same background shape 
parameter for the kinematic threshold function. Figure [9] 
shows an example fit for the combined B — > K*£'^£^ 
modes in the low s region. Table IVII summarizes the 
results. Figure [TO] shows Act as a function of s. Our re- 
sults are consistent with the SM expectation of negligible 
direct Acp- 

Wc fit the e+e~ and I^-^IJ.^ samples in the four K£'^£^ 
modes and four K*£^£~ modes in the low s and high 
s regions separately to measure the lepton-flavor ratios. 
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5.22 5.24 5.26 5.28 
m^g (GeV/c^) 



FIG. 5: The mEs spectrum in bin S4 for all Kt^l~ modes 
combined showing data (points with error bars), the total 
fit (blue solid line), signal component (black short-dashed 
line), combinatorial background (magenta long-dashed line), 
hadrons misidentified as muons (green dash-dotted line), and 
the sum of cross-feed and peaking components (red dotted 
line). 




5.26 5.28, 
m^g (GeV/c^) 



FIG. 6: The (a) tmes and (b) mKn mass spectra in bin si 
for all four K*t^t~ modes combined showing data (points 
with error bars), the total fit (blue solid lines), signal com- 
ponent (black short-dashed lines), combinatorial background 
(magenta long-dashed lines), hadrons misidentified as muons 
(green dash-dotted lines), and the sum of cross-feed and peak- 
ing components (red dotted lines) . 



TABLE V: Measured branching fractions [10 
s bin. The first and second uncertainties are 
systematic, respectively. 



^] by mode and 
statistical and 



s ( GeVVc'*) 



0.10-2.00 

2.00-4.30 

4.30-8.12 

10.11-12.89 

14.21-16.00 

>16.00 



1.00-6.00 



Slg 



B K£+i- 
B\W- 



20 

17.4 

37.1 



+5.9 
-5.4 
+5.4 
4.8 
+8.0 
7.5 



se.ol^ g 

19.7+_re 
22 s^^''^ 



0.71 
0.49 
0.94 
0.90 
0.49 
0.67 



+0.20 

0.18 
+0.15 
0.13 
+0.20 
0.19 
+0.20 
0.19 
+0.15 
0.14 
+0.23 
-0.21 



±0.02 
±0.01 

±0.02 
±0.04 
±0.02 
±0.05 



39.4 



+7.7 
7.1 



1.36 



+0.27 



±0.03 



B 



26.0 
14.5 
29.3 
31.6 
24.1 
14.1 



+ 7.1 
6.4 
5.3 
4.6 
+9.1 
8.3 
+8.8 
1 

+6.7 
6.0 
+6.6 
9 



1.89to 
0.95tg 
1.82^0 
1.86^0 



1.46"' 



1.02 



±0.06 
±0.04 
±0.09 
±0.10 
±0.06 
±0.06 



33.1 



h8.6 
-7.8 



2.05 



+0.; 
0.. 



Figure [TT] shows an example fit for the combined K^i'^ fi 
and Ke'^e^ modes in the high s region. Table IVIII and 
Fig. [Ushow Uk and Uk- for s > O.lGeVVc*. Our 
results are consistent with unity as expected in the SM. 

We fit the data in each s bin separately to determine 
Ai for the four combined K£'^£^ and four combined 
K*£~^i~ modes. Figure [T3] shows an example fit for bin 
32- The results are summarized in Table IVIIll Figure [T4l 
shows our measurements as a function of s in compari- 



> 
(3 



CO 

■o 



> 

0) 

O 



in 
OQ 



BABAR, 471 M BB 
CDF, 6.8 tb' 
Belle, 657 M BB 




2Q , 

s (GeV^/c^) 



FIG. 7: Partial branching fractions for the (a) K£'^£~ and (b) 
K*£^£~ modes as a function of s showing BABAR measure- 
ments (red triangles). Belle measurements ^2^] (open squares), 
CDF measurements [23| (blue solid squares), and the SM pre- 
diction from the All et al. model with B -)• if'*' form 
factors [U (magenta dashed lines). The magenta solid lines 
show the theory uncertainties. The vertical yellow shaded 
bands show the vetoed s regions around the J/ip and il){2S). 



K tr 





•B- 







BABAR, 471 M BB 
-■- CDF, 6.8 fb ' 

Belle, 657 M BB 
I I All '02 
I I Zhong 02 



x10"° 



0.5 



1.5 



Branching Fraction 



2.5 



±0.07 



FIG. 8: Total branching fractions for the A7+£" and K''-£+£- 
modes (red triangles) compared with Belle [2^ (open squares) 
and CDF [29I ] (blue solid squares) measurements and with pre- 
dictions from the Ali et al. [31 (light grey bands) , and Zhong et 
al. (dark grey bands) models. 
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son with those of Belle [28j. The two sets of results are 
seen to agree within the uncertainties. Our results are 
also consistent with the SM prediction that Ai is slightly 



negative {'^ —1%) except in bin si, where it is predicted 
to have a value around +5% jj]. 

Our Ai measurements in the low s region (0.10 < s < 
8.12 GeVVc") yield 




5.22 5.24 5.26 5.28, 
m^g (GeV/c^) 



0.9 1 1.1 

m^„(GeV/c^) 



FIG. 9: (a)&(c) rriES and (b)&(d) niKw fits for Acp in 
the (a)&(b) B and (c)&(d) B low s region for all four 
K*i^£~ modes combined. Data (points with error bars) are 
shown together with total fit (blue solid lines), combinatorial 
background (magenta long-dashed lines), signal (black short- 
dashed lines), hadronic background (green dash-dotted lines), 
and the sum of cross- feed and peaking background (red dotted 
lines). 



0.5 
0.4 
0.3 
0.2 
0.1 

-0.1 
-0.2 
-0.3 
-0.4 
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^ Kll' 
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FIG. 10: CP asymmetries Acp for Ki'^i' modes (red solid 
triangles) and K*£^£~ modes (red open circles) as a function 
of s. The vertical yellow shaded bands show the vetoed s 
regions around the J/tp and ^l)(2S). 



A^r{B 
A^r{B - 



K*£+£-) = 



-0.58t[5;3^ ± 0.02 [2.1cr], 
-0.25t[!:?? ± 0.03 [1.2cr], 



where the first uncertainty is statistical and the sec- 
ond is systematic. The Ai significances shown in the 
square brackets include all systematic uncertainties. We 
estimate the significance by refitting the data with Ai 
fixed to zero and compute the change in log likelihood 
v2Aln£ between the nominal fit and the null hypothe- 
sis fit. 



TABLE VI: Measured Acp by mode and s region. The 
first and second uncertainties are statistical and system- 
atic, respectively. "All" refers to the union of 0.10 < s < 
8.12 GeVVc* and s > 10.11 GeV^c*. 



s ( GcVVc*) 


Acp{B+ K+£+£-) 


Acp{B K*£+l-) 


All 


-0.03 ±0.14 ±0.01 


0.03 ±0.13 ±0.01 


0.10-8.12 
>10.11 


0.02 ±0.18 ±0.01 
-O.Oeloi? ±0.01 


-0.13io;w ±0.01 

o.ielo.lg ±0.01 




5.26 5.28, 
tTiEs (GeV/c^) 



5.26 5.28, 
iHes (GeV/c^) 



FIG. 11: rriES fits for Hk in the (a) Ke+e~ and (b) Kfj.+ ij,~ 
modes in the high s region. Data (points with error bars) are 
shown together with total fit (blue solid lines) , combinatorial 
background (magenta long-dashed lines), signal (black short- 
dashed lines), hadronic background (green dash-dotted lines), 
and the sum of cross-feed and peaking background (red dotted 
lines) . 



TABLE VII: Measured 7?.^(,) by mode and s region. The 
first and second uncertainties are statistical and system- 
atic, respectively. "All" refers to the union of 0.10 < s < 
8.12 GeV^c* and s > 10.11 GeV^c*. 



s {GeVyc") 


Hk 


Rk- 


All 


l-OOiois ± 0.07 


1.13lo;26 ±0.10 


0.10-8.12 
>10.11 


O.TAto.ti ± 0-06 
1.43t»ji±0.12 


1.06lo;33 ± 0-08 
1.18tJ! j^ ±0.11 
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X. CONCLUSION 

In summary, we have measured total and partial 
branching fractions, direct CP asymmetries, lepton-flavor 
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FIG. 12: Lepton flavor ratios TZj^{,) for the Kl^l~ (red solid 
triangles) and K*£^£~ modes (red open circles) as a function 
of s. The vertical yellow shaded bands show the vetoed s 
regions around the J/ip and ip{2S). 



ratios, and isospin asymmetries in the rare decays B — 
using 471 million BB pairs. These results 
provide an update to our previous measurements on 
branching fractions and rate asymmetries excluding the 
s < 0.1 GeV^/c^ region [32,]. The total branching frac- 
tions, B{B KI+1-) = (4.7 ± 0.6 ± 0.2) x 10"^ and 
B{B -> K*£+£-) = (10.2t};;^ ± 0.5) x 10~^ are mea- 
sured with precisions of 13% and 14%, respectively. The 
partial branching fractions as a function of s agree well 
with the SM prediction. For 0.10 < s < 8.12 GeV^/c^, 
our partial branching fraction results also allow compar- 
isons with SCET based predictions. CP asymmetries for 
both B K£+£- and B K*£+£- are consistent with 
zero and the lepton-flavor ratios are consistent with one, 
both as expected in the SM. The isospin asymmetries at 



TABLE VIII: Measured Ai by mode and s bin. The first 
and second uncertainties are statistical and systematic, re- 
spectively. 



s (GeVVc^) 



0.10-2.00 

2.00-4.30 

4.30-8.12 

10.11-12.89 

14.21-16.00 

>16.00 



1.00-6.00 



Ai 



B K£+£' 



B K''£+£'- 



-0.73 



0.95 
0.48 
0.55 
+0.27 
-0.30 
-+0.25 



±0.03 



omtoli ± 0.04 



-0.05Io 29 ± 0.03 



0.05 
-0.93 



-0.48lo;ig ±0.05 

-+0.61 



-4.99 



±0.04 



0.24+^;°^ ± 0.04 
l-07t^;9l±0.35 



-0.41 ±0.25 ±0.01 



-Q.2Qt°of3 ± 0.03 




5 10 15 2g 

FIG. 13: The rriES and mK-n fit projections for the (a)&(b) * (GeV I 

K*+£+£- and (c)&(d) K*'^£+£- modes in bin sa. Data 



(points with error bars) are shown together with total fit (blue FIG. 14: Isospin asymmetry Ai for the (a) K£^ £~ and (b) 

solid lines), combinatorial background (magenta long-dashed K*£^£~ modes as a function of s (red triangles), in compar- 

lines), signal (black short-dashed lines), hadronic background ison to results from Belle [28l | (open squares). The vertical 

(green dash-dotted lines), and the sum of cross- feed and peak- yellow shaded bands show the vetoed s regions around the 

ing background (red dotted lines). J /tj} and ^{2S). 
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low s values are negative. For 0.10 < s < 8.12 GcV^/c* 
we measure Ai{B K£+i-) = -0.58^° ± 0.02 and 
AiiB K*e+i-) = -0.25tg i? ± 0.03. The isospin 
asymmetries arc all consistent with the SM predictions. 
All results are in good agreement with those of the Belle, 
CDF, and LHCb experiments. 
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